Engineering Negative Differential Conductance with the Cu(lll) Surface State 
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Low-temperature scanning tunneling microscopy and spectroscopy are employed to investigate 
electron tunneling from a Ceo-terminated tip into a Cu(lll) surface. Tunneling between a Ceo orbital 
and the Shockley surface states of copper is shown to produce negative differential conductance 
(NDC) contrary to conventional expectations. NDC can be tuned through barrier thickness or Ceo 
orientation up to complete extinction. The orientation dependence of NDC is a result of a symmetry 
matching between the molecular tip and the surface states. 
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Important advances have been made over the past 
decades in molecular electronics. Among these is the 
demonstration that molecules can perform controllable 
functions such as negative differential resistance or con- 
ductance (NDC) [lH3[. First discovered in the Esaki 
diode Q , NDC leads to regions in the I ~V curve where 
the current / decreases (increases) with increasing (de- 
creasing) voltage V. This fundamental property is nowa- 
days exploited in CMOS devices for low-power memory, 
fast switches or oscillators. Early atomic-scale observa- 
tions of NDC by scanning tunneling microscopy (STM) 
have been attributed to narrow energy states in tip and 
sample analogy with the resonant tunneling lead- 

ing to NDC in semiconductors Q- In principle, resonant 
tunneling via molecular orbitals also leads to NDC in 
single molecules , but progress in this direction has 

been hindered by the lack of microscopic control over 
electrode and molecule status. 

It was shown in particular with single Cgo that NDC 
occurrence can be improved by narrowing the molecu- 
lar levels through a reduced molecule-substrate coupling 
111 Il2| . in order to exploit the bias-dependence of the 
transmission function l^, 14 1. In this Letter, we intro- 
duce a different approach and report the occurrence of 
single-molecule NDC with a Ceo-tcrminated tip. By at- 
taching a molecule to the STM tip an increased control 
is gained over the entire tunnel junction 15, [l^. This 
method allows, unlike previous NDC studies, exploring 
the NDC occurrence with well-defined pristine metal sur- 
faces serving as a conterelectrode. Taking advantage of 
this setup, we demonstrate in stark contrast to conven- 
tional mechanisms that NDC can be produced by elec- 
tron tunneling between a molecular orbital of the tip and 
a two-dimensional electron gas hosted by a copper sur- 
face — the Shockley surface states of Cu(lll). In this cal- 
ibrated setup, NDC may be tuned by varying the barrier 
thickness or by changing the Cgo orientation up to com- 
plete extinction. Our study demonstrates that molecular 
orbitals act as angular momentum filters to the tunnel- 
ing process, leading, in particular, to NDC if accurately 
matched with the local orbital symmetry of the surface 



states. These findings should simplify NDC engineering 
at the atomic-scale. 

The measurements were carried out with an STM op- 
erating at 4.6 K in the 10~^^ mbar range. After cleaning 
Cu(lll) by sputter/ anneal cycles, 0.15 monolayers of Ceo 
were dosed onto the surface from a crucible containing a 
99.9% pure powder heated to 400°C. Single Au atoms 
were instead evaporated onto the cold surface by heat- 
ing a high-purity gold wire. The differential conductance 
spectra. dI/dV(V) were acquired via lock-in detection 
with a bias modulation ranging from 3 to 10 mV rms 
at 670 Hz (sample bias is measured with respect to the 
tip). The etched W tip employed was cleaned by sput- 
ter/anneal cycles and treated in vacuo by soft indenta- 
tions into the clean surface. The tip apex was therefore 
coated with copper. 

Figure [H^a) presents the electronic structure of the 
Cu(lll) surface acquired with a metal tip (dashed line). 
Known spectroscopic features are evidenced [l3l- The 
steplike feature at —0.45 V is associated to Shockley sur- 
face states, which are an experimental realization of a 
nearly two-dimensional electron gas; the marked upturn 
below —0.9 V is due to copper d-states present below 
the (111) band gap. To engineer NDC with a Cgo tip, 
we attach a single Cgo to the copper-coated tip by re- 
peatedly bringing the tip into contact with the target 
molecule and pressing beyond contact [l^ . Figure [TJa) 
presents a conductance spectrum acquired with a Cgo tip 
over Cu(lll). The spectrum is nearly mirror symmetric 
to the spectrum acquired with a metallic tip positioned 
above a Ceo molecule [compare with Fig. S2(b) in the 
supplementary material]. The lowest unoccupied molec- 
ular orbitals are located at negative bias (LUMO: —0.5 V, 
LUMO+1: —1.9 V), while the highest occupied orbitals 
are at positive bias (HOMO: 2.0 V, not shown). The 
widths of the molecular states are typically 0.5 V and 
compare well with a Ceo molecule adsorbed on a metal 
surface 3"21|. Most importantly, NDC is detected be- 
low the LUMO+1. 

Among the numerous Ceo tips investigated, the ones 
producing NDC had a specific Ceo orientation relative 
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FIG. 1. (color online), (a)— (d) dl/dV over Cu(lll) for var- 
ious Ceo tips (1.0 V, 0.1 nA). The orientation of Ceo at each 
tip apex is determined by constant-current images over a Au 
atom (sample bias: —2.0 V, current: 0.1 nA, size: 12 x 12 A'^). 
Each orientation is flanked by a Hiickel simulation of the 
LUMO-l-1 orbital, and for better visibility, by the same or- 
bital cut along the xy plane (the LUMO-l-1 orbital is local- 
ized on the pentagons forming the Ceo skeleton). Panel (a) 
present also a dl/dV spectrum over Cu(lll) acquired with 
a copper-coated W tip (dashed line, feedback loop opened at 
-3.0 V and 5 nA). 



to the copper surface. Figures HJa) — (d) summarize our 
findings. In order to image the tip apex, constant-current 
images of a Au atom were acquired [15|. These images 
exhibit a molecular pattern matching the LUMO-(-l for 
biases close to —1.9 V. The orientations reported in Fig.[T] 
correspond to a Cgo adsorbed on a carbon atom [labeled 
corner in Figjlja)], on a pentagon-hexagon bond [5:6 in 
FiglTJb)], on a hexagon-hexagon bond [6:6 in FigHJc)] 
and on a hexagonal ring [hexagon in Fig|ljd)]. The apex 
Cgo was rotated by performing a tip contact with a Au 
atom [l3|. While the corner and the 5:6 adsorptions fa- 
vor NDC, other orientations tend to suppress it. We 
also found that the interface structure between the Cgo 
molecule and the copper tip has only a limited impact 
on NDC. Typical spectra acquired with different Cgo- 
tips, but with the same corner or 6:6 orientation are pre- 
sented in Figs, ma) and^b), respectively. Despite the 
Ceo-copper interface changes in these tips, as indicated 
by the shift of the LUMO and LUMO-hl (of the order 
of 0.5 and 0.2 V, respectively) [2l[, the dominant mech- 
anism leading to NDC is still the Cgo orientation. 

Along with molecular orientation, counterelectrode se- 
lection is extremely crucial to NDC. In our setup, the 
counterelectrode can be conveniently chosen by posi- 
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FIG. 2. (color online). dl/dV acquired above Cu(lll) with 
two different tips (1.0 V, 1.0 nA), having (a) same corner 
orientation for Ceo [as in Fig. [if a)], (b) same 6:6 orientation 
for Cao [as in Fig. (Tfc)]. 



tioning a Cgo tip above nanoscale objects supported by 
Cu(lll) or by changing the metal surface. Several coun- 
terelectrodes were inspected, either having a featureless 
electronic structure over the energy range of interest such 
as a Au atom on Cu(lll) or the Cu(lOO) surface, either 
presenting a marked electronic structure such as a cobalt 
islands on Cu(lll) or a Ceo molecule on Cu(lll). 
The dl/dV spectra, presented as supplementary mate- 
rial, never exhibited NDC. The fact that NDC is in- 
stead observed on Au(lll) as on Cu(lll) [2^, indicates 
that electron tunneling into the Shockley surface states 
is likely a key mechanism behind NDC. 

For completeness, we investigated the influence of 
barrier thickness on NDC in the line of recent studies 
1^ . I3I 24 1 . Figure |31[a) presents a typical set of spectra 
acquired at various distances, while Fig. quantifies 
the changes in NDC with distance (noted z) for a collec- 
tion of Ceo-tips. Through a tip displacement it is possible 
to tune NDC and, unlike previous work, even switch it 
off well within the tunneling regime at a distance of 3 A 
from surface. As established below, NDC extinction is a 
general property of our simplified setup not exclusive to 
Ceo- 

In the following we discuss the origin of the NDC with 
emphasis on the central role played by the Shockley sur- 
face states. To successfully engineer NDC, an accurate 
choice of the counterelectrode, of molecular orientation 
and of barrier thickness turns out to be crucial. To start, 
we adopt a WKB framework and focus on a junction 
comprising a sample metal surface and a Ceo tip with 
density of states (DOS) ps and pt, respectively. The zero- 
temperature current / as a function of sample voltage V 
is then 



I(y, z) cx 
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{E+V)pT{E)T{E,V,z)dE, (1) 



where T{E,V,z) = exp(-2^2(<I> - V/2 - E)z) is the 
transmission function in Hartree atomic units. The ex- 
plicit summation over magnetic quantum numbers m 
(projection of electron angular momentum on the com- 
mon symmetry axis) excludes the extremely small prob- 
ability for tunneling between states with different m 
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FIG. 3. (color online), (a) dl /dV spectra for various barrier 
thicknesses. The spectra are normalized by the conductance 
(G = I /V) at which the feedback loop was opened, respec- 
tively from top to bottom; 0.01, 0.1, 1 and 10 nS. (b) NDC 
intensity versus tip-sample distance (left vertical axis) and 
corresponding changes of G (right vertical axis). From the 
tunneling regime, we extract an apparent-barrier height of 
$ = 5 eV. The distance is estimated by taking z — for 
the contact conductance (G vs z curve acquired at 1 V), the 
dashed area highlighting the contact regime. The NDC inten- 
sity corresponds to the minimum of the dl/dV curves in (a). 
The lines correspond to the simulation detailed in the text: 
a realistic DOS (dashed line) and a stepped DOS without 
d-bulk contribution (dash-dotted line). 
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[lol 25 1 . For an accurate description of the tunneling 
current, we need just to distinguish between m = and 
m in the DOS. This distinction is necessary as, near 
r, the sp-like Shockley surface states of Cu(lll) only cou- 
ple to the TO = component of the tip states. To model 
the sample DOS, we thus take a stepped function cen- 
tered at —0.45 eV for the surface states and a low-energy 
exponential upturn for the copper bulk d-statcs to rep- 
resent [dashed line, negative half-plane. Fig. IH^a)] 
and a constant background plus exponential upturn for 
[dotted line, negative half-plane, Fig. ija)]. The 
total DOS (dash-dotted line) roughly mimics the elec- 
tronic structure evidenced in Fig. [ija). For the Cgo-tip 
we represent pt as a set of three Lorentzian peaks for the 
HOMO, LUMO, and LUMO-l-1 as described in the cap- 
tion to Fig. m^a). The m = (dashed line) and m 7^ 
(dotted line) components sum up (thick line) to ampli- 
tudes proportional to the free-molcculc degeneracies. 

As shown in Fig. lU^b), the corresponding simulated 
dl/dV captures the essential experimental facts. There 
is indeed NDC in the spectrum just below the LUMO-l-1 
down to distances close to z = 3 A. The role of the 
Shockley surface states in the NDC is evident from the 
energy-resolved current shown in Fig. |4l^c): At a sample 
voltage V = —2.35 V the surface state onset matches the 
LUMO-l-1 resonance [11,01. Upon a slight increase in the 
voltage the surface state step feature moves away from 
the molecular resonance. The associated loss in energy- 
resolved current (dark orange area) is only partially com- 
pensated by the corresponding shift of the sample chem- 
ical potential (light orange area), the tunnel current ef- 
fectively decreases thus giving rise to NDC [2^. The 
comparison to our experimental results is also satisfying 



FIG. 4. (color online), (a) DOS of low-symmetry Ceo-tip 
(upper half-plane, thick line) and Cu(lll) surface (lower half- 
plane, dash-dotted line: V = V, thick line: V = -2.35 V) 
used in the WKB simulation (m = and m 7^ components 
shown with dashed and dotted lines, respectively). Copper 
DOS: step function of unit amplitude centered at —0.45 eV 
in p™"" and a constant background set to 0.5 in p™^*^ with 
low-energy upturn. Ceo DOS: three Lorentzian peaks mim- 
icking the HOMO (-2.0 eV), LUMO (0.4 eV) and LUMO-M 
(1.9 eV) with widths 0.5 eV, 1.0 eV, and 0.5 eV, respectively, 
and amplitudes proportional to the free-molecule degeneracies 
(thick line). The occupied states are indicated by the shaded 
blue areas for the sample voltage V = —2.35 V where the 
Shockley surface state edge matches the LUMO-f 1 resonance, 
(b) Simulated dl/dV for (a) via Eq. ([T}. (c) Energy-resolved 
current for V — —2.35V (dark orange) and V — —2.40 V 
(light orange) for 2 = 6 A. (d) — (f) Similar to (a) — (c) but 
for a high-symmetry Ceo orientation where the m = com- 
ponent is suppressed for HOMO and LUMO-f 1. 



as shown by the dashed line in Fig. [Sfb), where we also 
include the NDC produced by a stepped function (dash- 
dotted line) without the d-bulk upturn. The increased 
NDC intensity without the d-states highlights that they 
are detrimental to NDC. The peculiar decay reported in 
Fig. O^b) is therefore the result of a competition between 
opposite effects for c?-states and Shockley surface states. 

Another salient experimental finding is the impact of 
Cgo orientation on NDC, which was set aside in the above 
discussion. In Fig. [1] it was found that NDC occurs with 
Ceo in low-symmetry configurations (e.g., corner orien- 
tation) but disappears in high-symmetry configurations 
(e.g., hexagon orientation). Our NDC originates from the 
Shockley surface states which select m = tip states. 
This suggests that for the high-symmetry Cgo orienta- 



4 





isolated 5:6/corner 6:6 


hexagon 


pentagon 




In 




C3V 


C51, 


HOMO 


Hu 


2A'+3A" A1+2A2+B1+B2 


A2+2E 


A2+E1+E2 


LUMO 


T\u 


2A'+A" A1+B1+B2 


Ai+E 


Ai+Ei 


LUMO+1 




A'+2A" A2+B1+B2 


A2+E 


A2+E1 



TABLE I. Coo orbitals in different local symmetries imposed by molecular orientation with respect to an isotropic surface. 



tions the m = component of the LUMO+1 vanishes. 
Indeed such an effect can be rationalized by considering 
the transformation properties of the Ceo tip orbitals us- 
ing group theory as follows. 

For isolated Ceo molecules of icosahedral symmetry Ih 
the HOMO is a fivefold degenerate state of Hu symme- 
try while the LUMO and LUMO-l-1 are triply degener- 
ate states of Tiu and Tig symmetry, respectively. When 
the molecule is adsorbed on a surface these degeneracies 
will be reduced and, importantly here, the orbital 
angular momentum components onto the common sym- 
metry axis perpendicular to the surface will depend on 
the molecular orientation. In Tab. U we list how the or- 
bitals are split in different relevant local symmetry en- 
vironments. For instance, a Ceo molecule adsorbed with 
a hexagon on an isotropic surface possesses a threefold 
rotation axis and can thus be considered as belonging to 
the C^v group. As a consequence it can be shown that 
the Tig LUMO-l-1 splits into A2+E representations [isj . 
However, from the character table of C^y it is evident 
that s, pz, and d^^ orbitals do not transform according 
to those representations. In other words, the LUMO-l-1 
does not have a m = component for the hexagon ori- 
entation. Contrarily, a Ceo molecule adsorbed on a 5:6 
bond or a corner atom only possesses bilateral symme- 
try as described by Cg, and the LUMO+1 is split into 
A'+2A". As s, Pz, and dz^ do transform in Cs accord- 
ing to those representations the LUMO+1 consequently 
contains a to = component in these low-symmetry con- 
figurations. 

More generally, based on Tab. U we can conclude that 
tunneling into the surface state is symmetry forbidden 
for orbitals which do not transform as the totally sym- 
metric representation (A'/Ai), i. e., the HOMO for pen- 
tagon and hexagons orientations and the LUMO+1 for 
6:6, pentagon, and hexagon orientations. Note also that 
there arc no such symmetry constraints for the LUMO or- 
bital. To simulate the dl/dV recorded with Ceo tips with 
high-symmetry orientation we therefore eliminated the 
component in HOMO and LUMO+1 [Fig. IDJd)]. 
As shown in Fig. |4l^e) — (f), indeed this leads to the 
disappearance of NDC as was observed experimentally 

[Fig.[irc)-(d)]. 

To summarize, we have shown that electron tunneling 
between a well calibrated molecule and a two-dimensional 
electron gas produces a tunable NDC. By controlling 



both tip and sample states and accounting for their sym- 
metries, we demonstrated the concept of local orbital 
symmetry matching [lo| . Barrier thickness or symmetry 
matching between the molecule and the counterelectrode 
can be used to switch the NDC on and off. As disper- 
sive two-dimensional electronic states are widely avail- 
able in metals and semiconductors, and moreover may 
be controUably modified through artificial nanostructures 
or molecular adlayers, our findings should simplify NDC 
engineering in atomic-scale metal-organic junctions. 

We thank G. SchuU, J. I. Pascual, W. A. Hofer, and 
D. Sanchez-Portal for fruitful discussions. 
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